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A  Study  of  Gas-Centered  Swirl  Coaxial 
Injectors  Using  X-ray  Radiography 


S.  Alexander  Schumaker1,  Alan  L.  Kastengren2, 


Stephen  A.  Danczyk1,  Malissa  D.A.  Lightfoot1, 

&  Chris  F.  Powell2 


•  Quantitative  spray  data  is  lacking  for  liquid-gas  rocket 
injectors,  especially  in  the  near-injector  region 

—  Quantitative  data  needed  to  develop  scaling  laws  and  validate  models 

—  Near-injector  region  is  of  critical  importance  since  this  is  typically  where 
combustion  occurs 

9  Most  liquid  rocket  injectors  operate  at  conditions  (high  flow 
rates  and  high  pressures)  which  produce  optically  dense 
sprays 

—  Laser  light  is  multiply  scattered  and  insufficient  useable  light  exits  the 
spray  to  allow  the  use  of  laser  diagnostics  in  the  near-injector  region 

—  Other  measurement  techniques  can  be  invasive,  qualitative  or  provide 
only  generalized  time-averaged  information 

•  Examples  include  mechanical  patternation,  ballistic  imaging  and 
traditional  shadowgraphy 

9  X-ray  radiography  is  currently  the  only  technique  that  allows 
quantitative  results  for  mass  distribution  in  the  near-injector 
region  and  inside  of  the  injector  itself 
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Other  Techniques 


*  Ballistic  Imaging  can  provide 
quantitative  information  on  the 
size  and  shape  of  intact 
structures  but  ignores  the 

droplet  field 


PDPA  can  be  used  on  these 
sprays  but  typically  not  in  the 
near  injector  region 


8HSA-90 


100 

Mechanical  Pattenization  is  the 
only  technique  that  can  directly  |io 
measure  mass  distribution  but  | 
again  not  in  the  near  injector  1 1 
region 

^  o.i 
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Previous  X-ray  Radiography  Studies 


9  To  date  most  of  the  x-ray  radiography  studies  have  focused  on 
diesel  and  gasoline  injectors 

—  Leick,  P.,  Kastengren,  A.L.,  Liu,  Z.,  Wang,  J.,  &  Powell,  C.F.,  11th  Triennial 
International  Conference  on  Liquid  Atomization  and  Spray  Systems,  Vail, 
Colorado,  July  2009 

—  Kastengren,  A.  L.,  Powell,  C.  F.,  Liu,  Z.,  Moon,  S.,  Gao,  J.,  Zhang,  X.,  and  Wang, 
J.,  22nd  Annual  Conference  on  Liquid  Atomization  and  Spray  Systems, 

Cincinnati,  Ohio,  May  2010 

9  Other  studies  include: 

—  Aerated-liquid  jets 

•  Lin,  K.C.,  Carter,  C.,  Smith,  S.,  and  Kastengren,  A.,  50th  At AA  Aerospace 
Sciences  Meeting,  Nashville,  Tennessee,  January  2012. 

—  Impinging  jet  injectors 

•  Halls,  B.  R.,  Heindel,  T.J.,  Meyer,  T.R.,  and  Kastengren,  A.L.,  50th  A I AA 
Aerospace  Sciences  Meeting,  Nashville,  Tennessee,  January  2012. 

—  Swirl-coaxial  injector 

•  Eberhart,  C.  J.,  Lineberry,  D.  M.,  Frederick,  R.  A.,  and  Kastengren,  A.  L., 

48th  AIAA  Joint  Propulsion  Conference,  Atlanta,  Georgia,  August  2012. 

9  First  time  used  in  a  spray  with  a  strong  gas  phase 
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Gas-Centered  Swirl  Coaxial  Injectors 


•  Why  Gas-Centered  Swirl-Coaxial  (GCSC)  injectors? 

—  Interest  exists  around  the  world  in  designing  new  LOX-Hydrocarbon 
rocket  engines  using  an  OX-rich  staged-combustion  cycle. 

•  Recent  papers  from  Russia,  China,  South  Korea  and  multiple  US  groups 

—  Successfully  used  in  Russian  rocket  engines  employing  a  staged- 
combustion  cycle  (e.g.,  RD-170) 

—  Such  a  design  and  typical  mixture  ratios  result  in  a  high  speed  gas  flow 
which  is  available  to  atomize  the  liquid  hydrocarbon  fuel 

9  Design  criteria  and  scaling  laws  are  needed  over  a  range  of 
conditions  and  geometries 

9  Previous  work  by  the  authors  using  cold  flow  techniques  have 
investigated  the  effect  of  numerous  geometric  variations  on  the 
size,  shape,  and  frequency  of  the  wall  bounded  film  and  spray 

—  Geometric  variations  include  cup  radius,  step  thickness,  gas  radius  at 
end  of  sheltered  lip,  initial  film  thickness,  cup  length  and  swirl 

—  Injector-wall  and  injector-injector  interactions  have  also  been  explored 

9  Cold  flow  techniques  include  laser  sheet  imaging  of  the  film, 
shadowgraphy  of  the  spray,  limited  PDPA 
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Gas-Centered  Swirl  Coaxial  Injector 


*  Tangential  liquid  inlets  create  a  swirling,  annular  fluid  which  is 
atomized  by  high-velocity  unswirled  gas 

9  Momentum  flux  ratio  (the  main  scaling  parameter)  was  varied 
between  60  and  145  with  4  different  geometries 

—  Geometries  represent  changing  the  number  of  inlets,  the  inlet  size  and 
the  inlet  area  (swirl  number) 

—  Outlet  size  also  examined,  but  not  presented  today 


4  tangential  inlets 


8  tangential  inlets 


rg  =rp=6.35  mm 
Lc=33.0  mm 
r0=  9.53  mm 
S=  1.52  mm 
t=1.65  mm 
djn=  1.6mm,  4H 
=1.6mm,  8HDA 
=0. 989mm,  8HSA 
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Nondimensional  Parameters 


The  relevant  nondimensional  parameters 
can  be  obtained  by  nondimensionalizing 
the  force  balance  for  a  surface  disturbance 
by  the  liquid  inertia  (P/V^r2)  and  multiplying 
by  T2/Adjst  to  eliminate  disturbance  parameters 
from  the  aerodynamic  force  term 


tlift 

Drag 


Assuming  worst  case  scenario  for  disturbance 
size  and  that  all  constants  are  of  order  1 : 


Centripetal 


Aerodynamic 
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Example  with  the  8H-ONPNTN  Geometry: 


Surface  Tension 

J_ 

We 


Force 

Aerodynamic 

®  total 

Liquid  Phase 
Viscous 

Centripetal 

Frc 

Surface 

Tension 

Value 

114 

0.1 

0.2 

0.02 

•  GCSC  injectors  operate  in  a  regime  where  aerodynamic  forces 
are  dominant 
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Momentum  Flux  Ratio 


O 


Total 


pA 

pA 


Name 

r0 

(mm) 

X 

(mm) 

rP 

(mm) 

s 

(mm) 

ODHNTN 

7.62 

1.65 

4.45 

1.52 

ONPNTN 

9.53 

1.65 

6.35 

1.52 

OUHUTD 

11.4 

1.32 

7.24 

2.87 

Tiotal 

*  Swirl  alters  the  shear  between  the  liquid  and  gas  phase  which 
must  be  taken  into  account  in  <J> 


*  Difficult  to  define  due  to  complex  compressible  flow  field 

*  Use  of  single  velocity  to  describe  flow  field  questionable 

*  Best  collapse  of  film  lengths  found  using  total  liquid  velocity 
and  a  calculated  gas  phase  velocity  and  static  density 

*  Total  liquid  velocity  calculated  using  measured  mass  flow  rate 
and  assuming  conservation  of  momentum  between  the 
tangential  liquid  injector  holes  and  liquid  cup  exit 

*  Using  the  measured  mass  flow  rate,  total  temperature  and 
static  pressure  in  the  gas  plenum  static  density  was  calculated 
assuming  ID  flow  and  a  calorically  perfect  gas 
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Flow  Laboratory 


*  Self  contained  mobile  system  capable  of  delivering  up  to  one 
kg/s  of  H20  &  GN2  at  pressures  in  excess  of  200  atm 

—  Requires  only  power,  LN2,  and  exhaust  from  host  facility. 

—  System  fully  rated  to  408  atm  (Allows  more  GN2  storage) 

—  Dedicated  Allen-Bradley  control  &  Pacific  Instruments  data  acquisition 
systems 

—  Fully  remote  operation 

—  High  speed  abort  system  on  all  data  channels  for  added  pressure  safety 
—  System  is  on  wheels  and  can  be  assembled  in  under  2  days 
—  Ran  almost  continuously  (24  hours/day)  for  two  weeks 
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X-Ray  Radiography  Details 


*  Experiments  performed  at  the  7-BM  beamline  of  the  Advanced 
Photon  Source 

9  Synchrotron  bending  magnet  to  produce  polychromatic  (white 
beam)  which  is  made  into  a  monochromatic  beam 

—  Monochromator  is  tunable  between  5.1  keV  to  12  keV 

—  10  keV  used  for  the  present  study  based  on  the  absorption  ot  the 
anticipated  test  conditions 

—  Beam  is  small  and  has  a  FWHM  of  5  x  6  microns 

—  X-ray  flux  of  2.5x1 011  photons/s/mm2 

9  Both  time-averaged  and  time-resolved  measurements  can  be 
made 

—  The  system  was  not  optimized  for  time-resolved  measurements 
as  only  time-averaged  ones  were  originally  planned. 

9  Beer’s  law  is  used  to  convert  these  measurements  into  an 
equivalent  path  length  of  water  (EPL) 

—  Path  length  is  the  average  across  each  measurement,  4  seconds  for 
time  averaged  results  and  1  microsecond  for  time-resolved. 
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Test  Setup  &  Matrix 


•  Normalization  by  the  the  intensity  of  incident  light  performed  in 
two  steps 


—  Titanium  foil:  beam  intensity  variations  during  a  scan 

—  Zero  absorption  case:  average  signal  level  from  the  5  points  in  the  scan  with 
the  highest  transmission  (Outside  the  spray-accounts  for  gas  absorption) 


Near  Injector  Optics  Layout 

Lxhaust 


Focusing  Mirrors 


*■ 


Titanium  Foil 


+> 


Pin  Diode 


Condition 

mg 

(g/s) 

m, 

(g/s) 

^total 

ra 

Res 

Re, 

We 

Frc 

4H-60 

46 

37 

58 

0.299 

380,000 

2,100 

42  0.158 

4H-90 

46 

30 

86 

0.299 

380,000 

1,800 

29  0.158 

4H-120 

46 

26 

117 

0.299 

380,000 

1,500 

21 

0.158 

4H-145 

46 

23 

141 

0.299 

380,000 

1,400 

17  0.158 

8HSA-60 

46 

32 

58 

0.261 

380,000 

2,200 

42 

0.162 

8HSA-90 

46 

26 

89 

0.261 

380,000 

1,800 

28  0.162 

8HSA-110 

54 

26 

107 

0.261 

440,000 

1,800 

28 

0.162 

8HSA-120 

46 

22 

120 

0.261 

380,000 

1,600 

20  0.162 

8HDA-60 

46 

49 

60 

0.406 

380,000 

2,100 

41 

0.145 

8HDA-90 

46 

40 

91 

0.406 

380,000 

1,700 

27  0.145 

8HDA-120 

46 

35 

118 

0.406 

370,000 

1,600 

21 

0.145 

8HDA-130 

33 

26 

128 

0.406 

260,000 

1,100 

12  0.145 

8HDA-145 

46 

31 

148 

0.406 

380,000 

1,364 

17  0.145 

Injector 
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Equivalent  Pathlength 


*  EPL  is  the  pathlength-integral  of  the  amount  of  water  in  the 
beam 

*  EPL  is  a  function  of  velocity,  droplet  size  and  mass  flux 
distribution 


70  r 


60  - 


«  50  - 


40  - 


E 

W  30 

O) 

1  20  - 


10  - 


Simulations  of  EPL  were  run  to  understand  the  effect  of  the 
velocity  and  mass  distribution  on  the  EPL  profile 

—  Simulation  assumes  axisymmetric  mass  &  velocity  profiles  &  divides  the 
spray  into  a  grid 

—  Droplet  size  is  assumed  uniform  across  the  spray 

—  Vertical  spacing  between  droplets  in  each  cell  is  set  to  achieve  correct 
mass  flux  and  velocity 


-20  -15  -10 


o 

r,mm 


10  15 


20 


8HSA-120 

(5mm  Downstream) 

Dd=75  pm 
rhL=35  g/s 
Vma=11.6rn/s 
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EPL,  V,  m/s 


EPL  Simulation 


r.mrn 


The  center  valley  in  EPL  can  not  be 
explained  by  a  mass  deficit  alone 

The  EPL  deficient  in  the  spray 
center  is  largely  caused  by  high 
center  line  gas  post  velocity 

—  300  m/s  gas  velocity  at  gas  post  exit 

—  Droplet  diameter  (Dd)  assumed  50  pm 
for  both  Cases 
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EPL,  p.m  EPL,  p.m 


Average  EPL  Profiles 


8HSA-120 
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Mass  Weighted  Velocity 


*  Droplets  across  the  entire  spray  are  undergoing  acceleration  in 
the  near  injector  region 

—  Acceleration  can  be  quantified  by  the  mass-weighted  velocity  (Vma) 

—  Vma  obtained  by  dividing  rh,  with  the  integral  of  the  average 
projected  liquid  density  profile 

—  Between  initial  gas-liquid  contact  and  3  mm  downstream 
acceleration  is  on  the  order  of  1 ,000  m/s2 

—  Between  3  mm  and  10  mm  acceleration  is  on  the  order  of  10,000 
m/s2 
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180°  Radial  Profiles 


*  Due  to  space  constraints  at  beamline,  testing  was  conducted 
with  the  spray  flowing  horizontal 

*  Two  scans  with  180°  rotation  to  investigate  potential  gravity 
effects 

—  180°  flipped  about  r=0  for  comparison 
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Momentum  Flux 


*  Uniformity  improves  with  increasing  momentum  flux  ratio  (<P)  up 
to  a  point 

—  Decrease  in  EPL  profile  as  O  increases  is  due  to  the  decrease  in 
the  liquid  mass  flow  rate 

—  Expected  from  earlier  work  /  views  downstream 


4H-ONPNTN,  5mm  downstream 
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EPL, 


Swirl 


0°  orientation,  5  mm  downstream 


Ra=0.299,  4H-120  RA=0.261 .8HSA-120  RA=0.406, 8HDA-120 
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Spray  Symmetry  &  rhL 


5mm  downstream 


mL=26  g/s 
4H,  Vgp=  295  m/s 
8HSA,  Vgp=  309  m/s 
8HDA,  Vgp=  250  m/s 
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Time  Resolved:  A  Single  Droplet 


Can  the  x-ray  radiography  provide  spray  statistics? 


—  Prior  radiography  results  have  not 
extracted  droplet  information 

Consider  the  simplest  case  of  a 
single  droplet  moving  thru  the  beam 

—  The  peak  of  the  path  length  corresponds 
to  the  droplet  diameter 


—  The  elapsed  time  of  departure  from  0  path  length  relates  to  the  velocity 


A  simple  approach  is  taken  to  get  droplet  size  and  velocity  from  the 
data  even  where  multiple  droplets  are  in  the  beam 


—  Each  peak  in  the  data  is  a  droplet 

—  Its  diameter  is  the  offset  from  the 
bounding  valleys  (minimum) 

—  Its  velocity  is  related  to  the  time 
between  bounding  valleys 

—  This  produces  known  biases,  but 
demonstrates  power  of  measurement 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


20 


Density 


Improving  Injector  Understanding 


*  Even  with  multiple  droplets  in  the  beam,  sizes  and  velocities  can 
be  extracted 


—  There  are  known  biases  in  the  technique  and  they  have  been  investigated 
and  an  improved  methods  is  in  the  works 

*  While  these  biases  mean  that  the  current  assessments  are  only 
semi-quantitative,  some  useful  results  exist 

—  Droplet  distributions  appear  to  be  log-normal  which  is  somewhat 
unexpected  and  likely  a  result  of  noise  level  in  the  data  and  biasing 

—  Changing  the  liquid  inlet  size  has  little  or  no  effect  on  the  atomization  and 
droplet-size  distribution 

—  However,  changing  the  liquid  inlet  area  (swirl)  has  a  large  effect  on  droplet- 
size  distribution 


Droplet  Diameter  (microns)  Droplet  Diameter  (microns)  Droplet  Diameter  (microns) 
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X-Ray  Conclusions 

^^^ra^aSiograpfiy^a^iseaTo  quantitatively  examine  tne  near 

injector  region  of  GCSC  injectors 

*  Time-averaged  measurements  quantified  a  significant  level  of 
asymmetry  not  previously  observed  in  these  sprays 

—  However,  the  mass  and  velocity  remain  coupled  and  interpretation  can  be  a 
challenge 

—  Indicates  that  engines  using  GCSC  injectors  must  be  able  to  tolerate  some 
level  of  asymmetry 

*  EPL  deficits  in  the  spray  center  were  shown  to  be  caused  by  high 
velocities  in  the  middle  of  the  spray 

*  Mass  averaged  velocities  show  high  acceleration  in  the  near 
injector  region  as  far  as  20  mm  downstream 

*  Spray  width  is  unaffected  by  d>  but  does  decrease  with  increasing 
swirl  level 

*  Time-resolved  radiography  might  be  able  to  help  decouple  mass  & 
velocity  by  providing  droplet  size  and  velocity  distributions 

—  A  simple  approach  demonstrates  the  utility  of  the  technique 

—  Improvements  have  been  identified  and  implemented  for  future  testing  to 
improve  fidelity 

*  X-ray  radiography  is  enabling  spray  insight  not  achievable  with 
other  available  diagnostics 
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Back  Ups 
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Time  Resolved:  A  Single  Droplet 


Can  the  x-ray  radiography  provide  spray  statistics? 


—  No  spray  statistics  possible  in  near  injector 
region  with  conventional  diagnostics  and 
these  dense  sprays 

—  Prior  radiography  results  have  not 
extracted  droplet  information 


Consider  the  simplest  case  of  a 
single  droplet  moving  thru  the  beam 


—  The  peak  of  the  path  length  corresponds  to  the  droplet  diameter 


—  The  elapsed  time  of  departure  from  0  path  length  relates  to  the  velocity 

*  A  simple  approach  is  taken  to  get  droplet  size  and  velocity  from  the 
data  even  where  multiple  droplets  are  in  the  beam 


—  Each  peak  in  the  data  is  a  droplet 

—  Its  diameter  is  the  offset  from  the  bounding  valleys  (minimum) 

—  Its  velocity  is  related  to  the  time  between  bounding  valleys 

—  This  produces  known  biases,  but  demonstrates  power  of  measurement 
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Droplet  Statistics 


Multiple  droplets  are  in  the  beam  at  any 
given  time 

Even  so,  droplet  sizes  and  velocities 
can  be  extracted 

—  There  are  known  biases  in  the  technique 
and  they  have  been  investigated 

—  An  improved  technique  using  curve-fitting 
is  in  the  works 

Droplet  distributions  currently  appear 
to  be  log-normal  which  is  unexpected 
and  likely  a  result  of  noise  level  in  the 
data  and  some  biasing 

Velocities  are  a  Log-Normal 
distribution,  but  this  is  more  typical 
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Path  Length/  Diameter  (|am) 


Comparison  to  Time-Averaged 


•  Time-averaged  results  can  be  related  to  droplet  diameters 
extracted  from  time-resolved  data 

—  Correlation  is  not  exact  because  periods  of  no  droplets  are  included  in 
the  time-averaging 

—  For  a  single  droplet,  the  mean  path  length  would  be  7iD/4 

*  Mean  diameter  is  LARGER  which  indicates  it  is  being 
overestimated  by  the  simple  procedure 


-  Could  be  a  result  of  upward 
diameter  bias 

-  Most  likely  the  result  of  not 
measuring  very  small  droplets 
(noise  and  resolution) 

-  Both  will  be  improved  upon  with 
the  new  DAQ  at  APS  and  new  data 
processing  here  at  AFRL 
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Velocity  as  a  Function  of  Diameter 


*  There  is  a  strong  correlation  between  droplet  diameter  and 
velocity  in  these  results 

*  At  least  some  of  this  is  artificial  and  a  result  of  the  current 
experimental  set-up 

—  Velocity  as  a  function  of  diameter  does  not  change  with  axial 
(downstream)  distance 

—  Larger  droplets  are  faster  and  have  a  wider  range  of  velocities,  which  is 
not  what  would  be  expected  for  a  shear-driven  spray 

—  Happens  because  droplets  must  be  in  beam  10  microseconds  to  be 
measured  as  an  independent  droplet  in  the  current  simple  assessment 


Mean  Droplet  Diameter  (microns) 
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Probability 


Improving  Injector  Understanding 


*  While  biases  mean  that  the  current  assessments  are 
only  semi-quantitative,  some  useful  results  exist 

*  Changing  the  liquid  inlet  size  has  little  or  no  effect 
on  the  atomization  and  droplet-size  distribution 

*  However,  changing  the  liquid  inlet  area  (swirl)  has  a 
large  effect  on  droplet-size  distribution 


0.08 


0.06 


0.04 


0.02 


.1 


ILarge  Inlets 
iSmall  Inlets 


lllllmiii 


Him 


100  200  300  400  500 

Droplet  Diameter  (microns) 


Droplet  Diameter  (microns) 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


29 


Simple  Procedure 


*  Automated  Matlab  process  was  developed  to  find  peaks  and 
troughs  and  translates  them  to  droplet  diameter  and  velocity 

9  Other  assumptions  include 

—  Droplets  are  spherical  and  travel  so  that  their  centers  are  captured  by 
the  beam 


—  Droplets  have  only  axial  velocities 

Signal-to-noise  ratio  not  optimized  here,  so  a  running-average 
is  taken  to  smooth  data 


—  Lower  x-ray  energy  will  be  used  in 
the  future  to  improve  the  signal-to- 
noise  and  hopefully  remove  the 
need  for  filtering 
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Obviously,  there  are  multiple  droplets  in  the  beam  at 
one  time  (troughs  are  not  0),  so  simple  procedure 
introduces  biases 

—  If  2  droplets  are  same  velocity  and  size  and  enter 
simultaneously,  the  procedure  finds  1  droplet  of  double  size 
and  velocity 

—  If  2  droplets  are  the  same  but  entry  is  offset  then  an  “extra” 
droplet  is  found  and  there  is  a  bias  to  underestimate 
diameter 

—  If  2  droplets  are  same  size  and  entry  time  but  different 
speeds  then  the  faster  droplet  appears  too  large  and  fast 
and  the  slower  droplet  has  an  overestimated  velocity 

Improvements  could  be  made  using  a  curve-fitting 
procedure  similar  to  in  spectroscopy 

—  However,  this  is  complex  and  time  consuming  (therefore,  it  is 
currently  incomplete) 


Distribution  A:  Approved  for  public  release;  distribution  unlimited. 


31 


Noise  Level 


*  The  automated  procedure  finds  “droplets”  in  the  case  where  no 
water  is  flowing 

—  Mean  diameter  is  23  microns  and  95th  percentile  diameter 
is  35  microns 


Droplets  are  also  found  outside  of  the  main  spray 

—  Spray  edge  is  difficult  to  determine  precisely 

—  These  values  are  larger  than  with  no  water,  95th  percentile  is  44-60 
microns  depending  on  the  test 


—  Mist  does  recirculate  in  the  exhaust  area;  this  can  be  mollified  by 
improved  exhaust 
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Probability  Mean  Diameter  (^m) 


Inlet  Size 


As  long  as  the  swirl  number  remains 
nearly  constant,  the  inlet  size  and 
number  should  not  effect  the  spray 

—  The  initial  shelter  should  allow  the  film  to  fully 
develop  and  have  no  “memory”  of  the  inlet 

A  slight  increase  in  mean  droplet 
diameter  is  observed  for  geometries  with 
larger  liquid  inlets 


—  Difference  is  small,  and  the  uncertainty  in  the 
results  is  not  fully  understood,  so  it  is  not 
clear  if  it  is  significant 

—  Swirl  is  very  slightly  different  between  the  two 
which  may  add  to  the  difference 

General  droplet  distributions  are  very 
similar  suggesting  similar  atomization 
behavior 
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Probability 


Swirl 


Shadowgraphy  shows  sprays  with  lower  swirl 
number  have  larger  droplets  on  the  periphery 

Mean  droplet  diameter  is  increased  at  lower  swirl 
number  (in  most  locations)  and  symmetry  is 
greatly  improved 

The  droplet  distribution  also  changes  suggesting  a 
difference  in  atomization  quality  and  behavior 
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